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The dynamical behavior caused by dry friction is stud- 
ied for a spring-block system pulled with constant velocity 
over a surface. The dynamical consequences of a general type 
of phenomenological friction law (stick-time dependent static 
friction, velocity dependent kinetic friction) are investigated. 
Three types of motion are possible: Stick-slip motion, con- 
tinuous sliding, and oscillations without sticking events. A 
rather complete discussion of local and global bifurcation sce- 
narios of these attractors and their unstable counterparts is 
present. 

PACS numbers: 03.20.+i, 46.30.Pa, 81.40.Pq 



Since more than 200 years Coulomb's laws of dry fric- 
tion have been well-known 0] . They state that the fric- 
tion force is given by a material parameter (friction coef- 
ficient) times the normal force. The coefficient of static 
friction (i.e., the force necessary to start sliding) is always 
equal to or larger than the coefficient of kinetic friction 
(i.e., the force necessary to keep sliding at a constant 
velocity) . 

The dynamical behavior of a mechanical system with 
dry friction is nonlinear because Coulomb's laws distin- 
guish between static friction and kinetic friction. If the 
kinetic friction coefficient is less than the static one stick- 
slip motion occurs where the sliding surfaces alternately 
switch between sticking and slipping in a more or less reg- 
ular fashion Q. This jerky motion leads to the everyday 
experience of squeaking doors and singing violins. 

Even though Coulomb's laws are simple and well- 
established (many calculations in engineering rely on 
these laws), they cannot be derived in a rigorous way 
because dry friction is a process which operates mostly 
far from equilibrium. It is therefore no surprise that de- 
viations from Coulomb's laws have often been found in 
experiments. Typical deviations are the following: (i) 
Static friction is not constant but increases with the stick- 
ing time (111, i- e -; tne ti me since the two sliding surfaces 
have been in contact without any relative motion, (ii) 
Kinetic friction depends on the sliding velocity; for very 
large velocities, it increases roughly linearly with the slid- 
ing velocity like in viscous friction. Coming from large 
velocities, the friction first decreases, goes through a min- 
imum, and then increases ||[| ■ In the case of boundary 
lubrication (i.e., a few monolayers of some lubricant are 
between the sliding surfaces) it decreases again for very 
low velocities (see Fig. |]) • The coefficient of kinetic 
friction as a function of the sliding velocity has therefore 



at least one extremum. The kinetic friction can exceed 
the static friction, but in the limit of zero sliding velocity 
it is still less than or equal to the static friction. 
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FIG. 1. Schematical sketches of typical velocity dependent 
kinetic friction laws for systems (a) without and (b) with 
boundary lubrication. 

The aim of this paper is to give a rather complete 
discussion of the nonlinear dynamics of a single degree 
of freedom for an arbitrary ■phenomenological dry fric- 
tion law in the sense mentioned above. This goes be- 
yond the discussion of specific laws found in the litera- 
ture |l|,||,§J§-§| ■ A phenomenological law for the friction 
force depends only on the macroscopic degrees of free- 
dom. This implies that all microscopic degrees of free- 
dom are much faster than the macroscopic ones. At the 
end of this paper I will give a simple argument why this 
assumption will not always be valid. To reveal this inva- 
lidity on the macroscopic level, it is therefore important 
to have a complete knowledge of the dynamical behav- 
ior under the assumption that this time-scale separation 
works. 

There are two other important reasons for knowing the 
consequences of the different dry friction laws: (i) Fric- 
tion coefficients can be measured only within an appara- 
tus (e.g. the surface force apparatus |n]] or the friction 
force microscope 0). Below we will see that the dynam- 
ical behavior of the whole system is strongly determined 
by the friction force and the properties of the apparatus. 
For example, stick-slip motion makes it difficult to obtain 
directly the coefficient of kinetic friction as a function of 
the sliding velocity. Thus, the influence of the measuring 
apparatus cannot be eliminated, (ii) Dry friction plays 
also an important role in granular materials p2[ . An 
open question there is whether or not the cooperative 
behavior of many interacting grains is significantly influ- 
enced by the dynamical behavior due to modifications of 
Coulomb's laws. 
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FIG. 2. A harmonic oscillator with dry friction. 

The mechanical environment (e.g. the apparatus) of 
two sliding surfaces may have many macroscopic degrees 
of freedom. The most important one is the lateral one. 
Here only systems are discussed which can be well de- 
scribed by this single degree of freedom. Figure || shows 
the apparatus schematically. It is described by a har- 
monic oscillator where a block (mass M) is connected 
via a spring (stiffness k) to a fixed support (see Fig. ||). 
The block is in contact with a surface which slides with 
constant velocity vq. The interaction between the block 
and the sliding surface is described by a sticking-time 
dependent static friction force Fs(t s tick) and a velocity 
dependent kinetic friction force Fk(v). For the equation 
of motion we have to distinguish whether the block sticks 
or slips. If it sticks its position x grows linearly in time 
until the force in the spring (i.e., kx) exceeds the static 
friction F S - Thus 



if \x\ < F s (t-t r )/K, 



(la) 



where t r < t is the time at which the block has stickcd 
again after a previous sliding state. If the block slips the 
equation of motion reads 

Mi + kx — sign (vq — x)Fk (|fo — (lb) 
if x ^ v or \x\ > F s {0)/k, 

where sign (x) denotes the sign of x. 

We start our investigation with Coulomb's laws of con- 
stant static and kinetic friction. As long as x < vq the 
system behaves like an undamped harmonic oscillator 
with the equilibrium position shifted by the amount of 
Fk/k. Thus there are infinitely many oscillatory solu- 
tions. Below we will see that some of them may survive 
in the case of velocity dependent kinetic friction. The 
equilibrium position of the block is x = Fk(vq)/ k. It is 
called the continuously sliding state. 

Every initial state which would lead to an oscillation 
with a velocity amplitude exceeding vo leads in a fi- 
nite time to stick-slip motion. Independent of the ini- 
tial condition the slips always start with x = Fs/k and 
x = vq. Thus the stick-slip motion defines an attrac- 
tive limit cycle in phase space. This is not in contra- 
diction with the fact that the system behaves other- 
wise like an undamped harmonic oscillator. The rea- 
son for that is that a finite bounded volume in phase 



space is contracted onto a line if it hits that part in 
phase space which is defined by ([la]). Stick-slip mo- 
tion requires a kinetic friction Fx which is strictly less 
than the static one. Usually the sticking time i S (i C fe = 
2(i*s — Fx) / (kvq) is much larger than the slipping time 
tsup = 2(tt - arctan[(F s - F K )v^ {kM)- 1 / 2 ]),JM/k'. 
The maximum amplitude of the stick-slip oscillation [i.e., 
maxt x(t)] is a monotonically increasing function of vq 
which starts at Fs/k for vq = 0. This is also true for a 
velocity-dependent kinetic friction force. 

The unmodified Coulomb's law lead to a coexistence 
of the continuously sliding state and stick-slip motion for 
any value of the sliding velocity vq. In the more general 
case of a velocity dependent kinetic friction this bistabil- 
ity still occurs but in a restricted range of vq. Especially 
there will be always a critical velocity v c above which 
stick-slip motion disappears. This is an everyday experi- 
ence: Squeaking of doors can be avoided by moving them 
faster. 

In order to be more quantitative we solve the equa- 
tion of motion for a linear dependence of Fx on v, i.e., 
Fk(v) — F KQ + 7U, with 7 > 0. Eq. ([fj]) becomes the 
equation of a damped harmonic oscillator which can be 
easily solved. Instead of a continuous family of oscilla- 
tory solutions we have an attractive continuously sliding 
state. Stick-slip motion disappears if the trajectory with 
x(Q) = Fs/k and x(0) = vq never sticks for t > 0. The 
critical velocity vo = v c is defined by x(t s u p ) — Fko/k 
and x(t s iip) = vq. It leads to two nonlinear algebraic 
equations for t s i ip and v c . For 7 <C \fM~K~ the solution 
can be given approximately: 



(2) 



The critical velocity v c plays an important role in the 
discussion of the nature of stick-slip motion, because 
its measurement tells us indirectly something about the 
mechanisms of dry friction (see the discussion in ||). 

Next we discuss a general nonmonotonic Fk(v) like 
the examples shown in Fig. [jj The static friction Fs is 
still assumed to be constant. The continuously sliding 
state exists for all values of vq but it is stable only if 
F' K = dFK{vo)/dv > 0. At an extremum of Fk(v) the 
stability changes and a Hopf bifurcation occurs. Near the 
extremum and for small deviations from the continuously 
sliding state the dynamics of 
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is governed by the amplitude equation (normal form) |T^ 1 
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If the third derivative of the kinetic friction at an ex- 
tremum is positive, the Hopf bifurcation is supercritical, 
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and in addition to the well-known attractors mentioned 
above, another type of attractor appears. I call it the 
oscillatory sliding state. It is a limit cycle where the 
maximum velocity remains always less than vq. Thus 
the block never sticks. Its frequency is roughly given by 



the harmonic oscillator of the left-hand side of (lb). The 
second derivative of the kinetic friction is responsible for 
nonlinear frequency detuning. Note that the frequency of 
the stick-slip oscillator is usually much smaller than the 
frequency of the oscillatory sliding state. This oscillatory 
state is similar to the limit cycle of Rayleigh's equation 
u + e(ii 3 — u) + u = |0|, in fact, Rayleigh's equation is 
a special case of (ffj|). Depending on the kinetic friction, 
several stable and unstable limit cycles may exist. By 
varying v they are created or destroyed in pairs due to 
saddle-node bifurcations. 

It should be noted that the Hopf bifurcation described 
by ([|) is not related to the Hopf bifurcation observed by 
Heslot et al. j| which occurs in a regime (called creeping 
regime) where (0) is not applicable (see also the discus- 
sion below about the validity of dry friction laws). 

An oscillatory sliding state exists only if its maximum 
velocity is smaller than the sliding velocity vq because of 
the sticking condition (p"aj). How does the interplay of 
the oscillatory sliding states and the sticking condition 
leads to stick-slip motion? In order to answer this ques- 
tion we calculate the backward trajectory of the point 
lim £ ^o(-F'fs"(0)/K, vo — e) in accordance with (|Tb|). Three 
qualitatively different backward trajectories are possible: 

1. The backward trajectory hits the sticking condi- 
tion. Together they define a bounded set of initial 
conditions leading to non-sticking trajectories. I 
call the boundary of this set the special stick-slip 
boundary; it is not a possible trajectory but it sep- 
arates between the basins of attraction of the stick- 
slip oscillator and the non-stick-slip attractors. 

2. The backward trajectory spirals inwards towards 
an unstable oscillatory or continuously sliding 
state. Again all initial states outside these repelling 
states are attracted by a stick-slip limit cycle. 

3. The backward trajectory spirals outward towards 
infinity, and stick-slip motion is impossible. 

Two types of local bifurcations arc possible: If the 
backward trajectory changes from case 1 to case 3 the 
stick-slip limit cycle annihilates with the special stick- 
slip boundary. For changes from case 1 to case 2 the 
special stick-slip boundary is either replaced by an un- 
stable continuous or oscillatory sliding state or it annihi- 
lates with a stable continuous or oscillatory sliding state. 
A change from case 2 to case 3 is not possible. Fig- 
ure [| shows for a particular choice of Fk{v) both types 
of bifurcations. Here the first bifurcation type occurs at 
v w 0.059, 0.082, and 0.966. The second type occurs at 



vq 0.162 and 0.785. This example shows that for in- 
creasing vq stick-slip motion can disappear and reappear 
again. 

Besides of the well-known bistability between stick-slip 
motion and continuous sliding j^] , multistability between 
one continuously sliding state, several oscillatory sliding 
states, and one stick-slip oscillator is possible (see Fig. |^). 
Eventually for large sliding velocities all attractors except 
that of the continuously sliding state will disappear be- 
cause the kinetic friction has to be an increasing function 
for sufficiently large sliding velocities. 
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FIG. 3. Typical bifurcation scenarios for a partic- 
ular kinetic friction force Fk(v) of type shown in 
Fig. ^(b). The following function has been chosen 

F K (v) = v[ji +72 + (72 - 7i)(« - v)/y/(Av) 2 + ( v ~ *) 2 ]A 
with 71 = 3, 72 = 0.1, v = 0.2, and Av = 0.05. The results 
are obtained by numerical integration of the equation of mo- 
tion (Q). The other parameters are Fs = 1, M = 40, and 
K = 1. Solid (dotted) lines indicate stable (unstable) contin- 
uously sliding states (CS), oscillatory sliding states (OS), or 
stick-slip motions (SS). The dashed-dotted line indicates the 
special stick-slip boundary. 

The strongly overdamped limit (i.e., \dFj( (v) / dv\ ^S> 
V kM for any v except in tiny intervals around the ex- 
trema) leads to a separation of time scales: From an 
arbitrary point [x, x) in phase space with x < vq the sys- 
tem moves very quickly into the point (x, v) where v is a 
solution of kx = Fk(vq — v) with F' k (vq — v) > 0. Points 
on the curve kx = Fk(vq — x) with F' K < are unstable. 
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They separate basins of attraction of different solutions 
v. After the fast motion has decayed the system moves 
slowly on the curve kx = Fk(vq — x). The direction is 
determined by the sign of x. It either reaches a stable 
continuously sliding state, or, near an extremum of Fk, it 
jumps suddenly to another branch of the curve or to the 
sticking condition. For kinetic friction laws of the form 
shown in Fig. |l|(b) with vq between the two extrema, we 
get an oscillatory sliding state. It is a relaxation oscilla- 
tion which may be difficult to distinguish from a stick-slip 
oscillation. In the case of a friction law with a single min- 
imum at v = v m as shown in Fig. |l](a) we get stick-slip 
motion for vo < v m Q . In the strongly overdamped limit 
any multistability disappears except near the extrema of 
Fk (v) . The experiments of Yoshizawa and Israelachvili 
fl4| are consistent with the assumption that the system 
is in a strongly overdamped limit with a friction law as 
shown in Fig. 0(a) @. 

In order to discuss the influence of a stick-time depen- 
dent static friction on the stick-slip behavior we define a 
stick-slip map x n +\ = T(x n ), where x n is the position of 
the block just before slipping. For constant static fric- 
tion the map reads T{x) = Fs/k. The position just at 
the time of the slip-to-stick transition is defined by x s n . 
It is a function of x n , i.e., x s n = g(x n ), where g is usually 
a monotonically decreasing function. The sticking time 
fstick j g smallest positive solution of 



F s {tT ck ) 
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This defines a function tf^ zck = h(x s n ) which is a mono- 
tonically decreasing function due to F' s > 0. Thus the 
stick-slip map is given by T{x) — Fs(h(g(x)))/ k. If the 
map has one fixed point, then stick-slip motion exists. 

For Fk = const = F$(0), stick-slip motion disappears 
if v a > v c — sup t>0 2[Fs{t) — Fs(0)]/{nt). For a non- 
convex Fs{t) the supremum occurs at a non-zero value 
of the sticking time leading to a saddle-node bifurcation 
of a stable and an unstable fixed point of the stick-slip 
map. At vo = v c the stick-slip motion has a finite ampli- 
tude, contrary to the case of a convex Fg(t) |8). Because 
T is a monotonically increasing function, limit cycles or 
even chaos are not possible. If the slip-to-stick transi- 
tion does not happen at the first time when x becomes 
equal to vo [because of |x*| > Fs(0)/k] chaotic motion 
may occur Q. In this case we get a non-monotonic T 
due to a non-monotonic g. Such over-shooting is only 
possible if Fs(oo)/Fs(0) becomes relatively large. For 
example, for a constant kinetic friction over-shooting oc- 
curs if F s (oo)/F s (0) > 1 + F K (0)/F s (0). For most re- 
alistic systems this condition is not satisfied. Note that 
the possibility of chaos is not in contradiction with the 
fact that the equation of motion ([[]) with constant F$ 
can not show chaotic motion. But the retardation of F$ 
turns ([fl) into a kind of differential-delay equation. 



Using phenomenological dry friction means that we 
treat dry friction as an element in a mechanical circuit 
with some nonlinear velocity-force characteristic like, say, 
a diode in an electrical circuit. This treatment is justified 
as long as the macroscopic time scales are much larger 
than any time scale of the internal degrees of freedom of 
the interacting solid surfaces. But there is one internal 
time scale which diverges if the relative velocity between 
the surfaces goes to zero: It is given by the ratio of a 
characteristic lateral length scale of the surface and the 
relative sliding velocity. Thus any kinetic friction law 
Fk{v) becomes invalid if 

^ microscopic length scale 

v ~ : : r~ ■ (6) 

macroscopic time scale 

The characteristic length scale ranges from several mi- 
crometers to several meters. It may be the size of the 
asperities, the size of the contact of the asperities, the 
correlation length of surface roughness, or an elastic cor- 
relation length. This limitation of any dry friction law 
does not concern oscillatory sliding states and continu- 
ously sliding states, as long as their relative sliding ve- 
locity stays always much larger than the critical velocity 
(^). But the transition between sticking and sliding in 
a stick-slip motion may be strongly affected by the fact 
that just after stick-to-slip transitions and just before 
slip-to-stick transitions, details of the interface dynamics 
become important. One may expect that the importance 
of these details increases when the maximum slipping 
velocity decreases. For example, Heslot et al. [|| found 
experimentally a completely different behavior when the 
maximum relative sliding velocity during a slip was below 
the critical value (§|). 

In this paper the nonlinear dynamics of a harmonic 
oscillator sliding over a solid surface has been discussed 
under the assumption that dry friction can be described 
by a velocity dependent kinetic friction and a sticking- 
time dependent static friction. Besides of the well-known 
continuously sliding state and the stick-slip oscillator, an 
oscillatory sliding state without sticking has been found. 
All typical bifurcation scenarios of these states are shown 
in Fig. | 
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